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Ion Charge States in the Fast Solar Wind: New Data Analysis 

and Theoretical Refinements 

J. Martin Laming^ & Susan T. Lepri^ 

O 

O ■ ABSTRACT 
(N 

D I We present a further investigation into the increased ionization observed in 

^ ' element charge states in the fast solar wind compared to its coronal hole source 

' regions. Once ions begin to be perpendicularly heated by ion cyclotron waves and 

execute large gyro-orbits, density gradients in the flow can excite lower hybrid 
waves that then damp by heating electrons in the parallel direction. We give 
further analysis of charge state data from polar coronal holes at solar minimum 



CN ■ and maximum, and also from equatorial coronal holes. We also consider further 

^ ' the damping of lower hybrid waves by ions and the effect of non-Maxwellian 

O ■ electron distribution functions on the degree of increased ionization, both of which 

appear to be neghgible for the solar wind case considered here. We also suggest 
that the density gradients required to heat electrons sufficiently to further ionize 
i-( • the solar wind can plausibly result from the turbulent cascade of MHD waves. 



X 



Subject headings: Sun: solar wind — atomic processes — plasmas — waves 



1. Introduction 

The mechanism(s) by which electrons and ions can exchange energy collisionlessly, that 
is in the absence of Coulomb collisional coupling, have become important issues in a variety 
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of astrophysical contexts. The fast solar wind provides an excellent laboratory to study such 
mechanisms in detail due to the relatively large amount data available. During its initial polar 
passes, in-situ observations from the Solar Wind Ion Composition Spectrometer (SWIGS) 
on Ulysses revealed that the fast solar wind is more highly ionized than would be expected 
given the electron temperature and charge states spectroscopically observed in its coronal 
hole source regions. Ulysses SWIGS observed charge states re lated to electro n temperatures 
near 1.5 x 10^ K, assuming collisional ionization equilibrium (IKo et al.lll997l ). contradicting 
the finding from SUMER on SOHO which indicate electro n temperatures near 8 x 10^ K out 
to a heliocentric distance of 1.5 Rq ( iWilhelm et al.l 119981 ). For Fe, tempera tures of 8 x 10 ^ 



K in the freeze-in regioii would result in a main charge state of 8+ (see e.g. iMazzotta et al. 



19981 : iBryans et al.ll2006l ). which is rarely observed in any substantial abundance in the solar 
wind. Gharge states in the fast solar wind for Fe are typically near 10+. Hence, any heating 
mechanism for the electrons between 1.5 solar radii and 3 solar radii needs to explain the 
discrepancy between these SUMER observations and the Ulysses SWIGS data. 

The fast solar wind is likely accelerated by io n cyclotron waves whi c h damp by heat - 



ing ions to high perp endicula r temp eratures (e.g. iGranmer et al.l Il999bl : iGranmer I l2000l ). 
In an initial analysis, iLaming] (120041 ) pursued the idea that some of this energy "leaks" to 
the electrons in the solar wind acceleration region by a collisionless process. The increased 
electron temperature then further ionizes the plasma. The most likely collisionless mecha- 
nism appeared to be the excitation of lower-hybrid waves by ions (most likely a-particles) 
executing large gyro-orbits in a cross- field density gradient. In this way, the electrons are 
not heated until the ions develop large gyror adii, beginning at abo ut 1.5 Rq heliocentric dis- 
tance. Thus the observational constraint of IWilhelm et al.l (119981 ) of electron temperatures 
of around 8 x 10^ K out to 1.5 Rq is not violated, which might not be the case for other 
means of electron heating. 

While the existence of such density g radients is hi ghly plausible from Interplanetary 
Scintillation observations (see discussion in iLamingl 120041 ). the origin of these structures was 
not addressed. The turbulent cascade of MHD Alfven waves excited at the solar surface is 
known to be essentially 2D in character, proceeding much faster in the direction perpen- 
dicular to the magnetic field than in the direction parallel to it. This means that direct 
damping of the Alfvenic turbulence on ions is very slow, because the Doppler shift, k\\v\\, is 
small and few ions are resonant. A possible solution to this lies in the inclusion of fast mode 



waves in the turbulence (IGhandran 



20051 : ILuo fc Melrosd l2006l ) . In such cases, the cascade 



is not confined to the perpendicular direction and high wave frequencies are accessible. Fast 
modes are also compressible, and more likely to produce the density gradients observed in 
Interplanetary Scintillation, and also invoked in thi s work to provide t he instability that gen- 
erates lower hybrid waves. In another approach, (IMarkovskiil l200ll : iMarkovskii &: HoUweg 
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20021 : iMarkovskii et al.ll2006l ) suggest that the perpendicular gradients of density, magnetic 
field or velocity produced by the turbulent cascade then generates cyclotron resonant waves 
by a plasma instability, which provide the minor ion heating. We argue that a similar plasma 
instability sh ould produce lower hybr id waves that damp on the electrons in the range be- 



yond 1.5 Rq. IMarkovskii et al.l (120061 ) consider in most detail the velocity shear produced at 
the high k± end of the turbulent cascade, this being the most effective means of generating 
ion cyclotron waves. Hence aside from being an interesting curiosity with interesting impli- 
cations in various fields of astrophysics, our interpretation of the increased ionization of the 
fast solar wind in terms of small scale perturbations, in our case density gradients, may have 
an important bearing on the central problem of the acceleration of the fast solar wind. 

In this paper we provide further analysis of in situ solar wind ionic composition data 
from Ulysses, for fast wind from polar coronal holes, and the Advanced Composition Explorer 
(ACE) for low latitude fast wind streams. These measurements generally support the earlier 
results, though the degree by which the ionization is enhanced appears to be lower than 
previously reported. We also describe some refinements to the theory, which combined with 
the lower observed enhancement in ionization, make it plausible that the cross field density 
gradient required is the end result of the MHD turbulent cascade. This pa per is or g anize d 



as follows. Sections 2 and 3 give some refinements to the theory presented in lLamingl (120041 ). 
Section 2 considers ion heating by lower hybrid waves, and show that this is always negligible 
compared to the electron damping of the wave for solar wind conditions. Section 3 discusses 
the electron distribution function arising from heating by lower hybrid waves, and shows both 
on theoretical and observational grounds that a Maxwellian is the most likely distribution. 
Section 4 presents new measurements of the Fe charge states in fast solar wind from polar 
coronal holes observed by Ulysses and equatorial coronal holes observed by ACE, section 5 
gives some further discussion and section 6 concludes. 



2. Ion Heating by Lower Hybrid Waves 



In the model of iLamingl (120041 ). lower hybrid waves are excited by a-particle execut- 
ing large gyro-orbits in the presence of a cross field density gradient. These waves are 
electrostatic ion oscillations with wavevectors concentrated in a narrow range (typically 
—ujpi/ujpe < COS 9 < ujpi/ujpe, where Upe and ujpi are the electron and ion plasma frequencies 
respectively) around the direction perpendicular to the magnetic field. Electron screening, 
that would normally inhibit the wave motion, is restricted for wavelengths greater than the 
electron gyroradius {k± < Ve/^e)- Since u!/k± >> uj/k\\, the wave can be simultaneously 
in Landau resonance with ions moving perpendicularly to the magnetic field, and electrons 
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moving parallel. Hence a wave generated by a minor ion executing a large gyro-orbit may 
in principle be damped by ions (most likely by protons) as well as by electrons. Ions may 
Landau damp the wave, but above a critical wave electric field, E = B {VLi/ojf^ oj/Ak^c 
(jKarneylll978l ). stochastic heating also sets in. We estimate the heating rates for ions and 
electrons as follows. Th e velocity diffusion coefficient for unmagnetized ions in electrostatic 
waves is (lMelrosdll986l . equation 6.43) 
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where R is the ratio of electric to total energy in the wave, ykj is the number density of 

wave quanta with wavevector k, e is the wave polarization vector and v is the ion velocity. 
Writing E^/8n = J RN (k) hwdP / (27r)^ we obtain 



-6 (tu — kv cosi/j) 



(2) 



Following iKarneyl (Il979l ). the diffusion coefficient for magnetized ions in perpendicularly 
propagating electrostatic waves is obtained from this result by averaging over Larmor angle 
■ip, with k — > k±, V — *• v±, 
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where we have assumed two resonances per gyro-orbit, i.e. the 5 function is satisfied twice 
in 27r. The heating rate is 
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where the second step follows from writing Of /dt = {l/v±) d/dv± {vj_Ddfi/dv±), and inte- 
grating by parts in cylindrical coordinates. Substituting the zero field ion velocity diffusion 
coefficient, the heating rate 
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where we have put the wave energy density W = uj\Eo\ /Sn x deL/dou = u^^ \Eq\ jAmuP'. 
The final step may be rewritten Q = 2WjLD/ni, where 'jj.n = y/7rtox^ exp —x"^ with x = 
ijjj ^/2k±Vt± is the ion Landau damping rate (see e.g. lLaminej|200ll . equation A12 with U = 0). 
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In nonzero magnetic field, we substitute equation 3 into equation 4, and with 

^i/2ft±, we get 
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with V = u^. The l ast form of the integral in equation 6 may be evaluated in terms of a 
Whittaker function JOradshtevn fc R.vziklnJ . 3.383.4, p319) to give 

Q = "^ux'T f 1) exp (-^^) Wy,,_y, {x^) . (7) 



This is then evaluated using lAbramovitz fc StegunI (119701 ) equation 13.1.33 (plQO) to give 

Q = —ujx^^exp -x^U ( -i, 0, j . (8) 



Using limiting forms for the confluent hypergeometric function U (—1/2, 0, x^) (lAbramovitz fc Stegun 



1970 



equations 13.5.11, pl93 and 13.5.2, pl93) we finally arrive at 
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i) ux^^/n exp —X 

Thus in the limit x >> 1, the r nagne tized ions are heated a factor of x/2 faster than the 
unmagnetized ions. iBrambillal (119981 ) gives a similar result (page 504), but his diffusion 
coefficient is different from our equation 3 by a factor Qi/2k±v±, where fli is the ion gyrofre- 
quency. For x ^ 0, the unmagnetized Landau damping is a factor 2^/^^ faster. 



The electron Landau damping rate is given by equation All in iLamingl (120011 ) as 
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which gives a damping rate of order A/vrcumj/mey^ exp — y2/2 where y = uj j ^2k\^V(,^. For 
X ~ ?/, this is mi/ me faster than ion Landau damping, and approximately ^/mi/m^ faster 
than the stochastic heating rate for magnetized ions. We conclude that lower hybrid waves 
damp almost exclusively by heating electr ons, and that the electric field in the wave may 
be higher than that set by iKarneyl (119781 ). For ion heating to be significant, we require 
y » X » 1, i.e. the ions should already be significantly hotter than the electrons, but 
still in a "cold plasma" regime. The only qualification to this might be if the lower hybrid 
waves propagate at cosines cos^ 6 << m^/mi so that k\\ —>■ 0. This seems unlikely, since even 
in the cold plasma case where lower hybrid growth rates and frequencies remain nonzero at 
cos^ = 0, they are still lower than at cos^ ^ = mg/mj, and so the existence of waves with 
nonzero kti and the consequent electron heating would seem to be inevitable. 
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3. Electron k Distributions in the Fast Solar Wind? 

Electrons heated by colhsionless processes do not necessarily maintain a Maxwellian 
velocity distribution function. In this section we investigate the effect of a nonMaxwellian 
electron velocity distribution function on the ionization balance of oxygen. The electrons 
are taken to be be heated into a k distribution, /(f) = {I + v"^ /2Kvf) As k ^ oo, 
/ (v) — > exp (— f ^/2fj^) and a Maxwellian is recovered. In Appendix B we give a simple 
derivation of this form for the electron velocity distribution function from lower hybrid wave 
heating. In the case that to « ^/2kvi, a k function results, with the value of n dependent 
on the wave electric field and the ambient density. In the opposite limit, lu >> \/2kvi, a 
Maxwellian results, with quasi-thermal velocity also dependent on the wave electric field. 
Based on our knowledge of coronal parameters, we would expect the latter case to be most 
likely for the acceleration region of the solar wind, where the ions are further ionized. In 
this section we show that the electrons in a Maxwellian distribution function are also more 
consistent with the observed charge states than are k distributions. 

We assume that the electron n distribution is isotropic, mainly for ease of computation. 
We have 

so = (1 — 5/2/t) fcT/m. The energy in the "tail" portion of the distribution function 
is (3/2) {nkT — nmv^) = IbnkT / An. We write the temperature equilibration rate between 
the "core" and the "tail" as 5T/2Krs where Ts = 27mv^ /SUp^lnA is the stopping time by 
Coulomb collisions for fast electrons in cold plasma. 

We modify the ionization and recombination rates according to the departure from the 
Maxwellian distribution function. In principle, the cross section for each process should 
be integrated over the k distribution to evaluate the new rate. This is straightforward 
for direct collisional ionization and radiative recombination, but for the resonant processes 
of dielectronic recombination and autoionization following inner shell excitation, becomes 
considerably more cumbersome. Instead, we approximate the k distribution by 

/ (.) = (1 + „72k.?) - = -P (^) + i exp + ... (12) 

which has been found by trial and error. We illustrate the comparison between the true and 
approximate k distributions in Figure [H We then evaluate approximate rate coefficients by 
making the appropriate sum over Maxwellian rate coefficients, and proceed with simulations 



- ^ i2nv^nf' ^ = InkT, (11) 
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in the manner described in iLamingl (120041 ) . We keep the same definition of the electron-ion 
equihbration parameter, 7' = jiMi/ujAfq'^ x {uj'^/k'^v±) where 7j is the lower hybrid wave 
growth rate due to the gyrating ions, Mj is the ion mass, A is the ion abundance, / the 
charge state fractions, q the charge, u the lower hybrid wave frequency, k the wavevector 
and v± the ion gyration velocity. Th e ion-ele c tron e nergy transfer rate is modified by a factor 
x2a;pg/f2g. This corrects an error in iLamind (120041 ) where only the electric energy density of 
the lower hybrid wave was considered in computing the energy transfer as 2H^7e, whereas W 
should properly be the total wave energy density. Thus f or a given 7^ t he electron internal 
energy, 3nkT/2, is the same in these models as in those of iLaming (120041 ). multiplied by this 
correction factor of 2ujp^/Ql ^ 20. Each simulation is initiated with a high value of k = 10^ 
to enforce an initial Maxwellian distribution close to the sun. After every time step in which 
the electron and ion temperatures and charge states are integrated, we update vf along with 
the other hydrodynamic variables in operator splitting fashion, according to 



dvl 
dt 



5 kT 

2k, mr^ 



(13) 



and then update k = 2.5/ (1 — mv^ /kT). This assumes that all of the electron heating goes 
to electrons in the tail. The Coulomb equilibration between the "tail' and the "core" is 
slower than the overall coUisionless heating rate, and so k decreases as the heating proceeds. 

Table 1 gives values of the 0+^/0+^ charge state ratio for models with varying 7', for 
initial wind fiow speeds of 5, 10, and 20 km s^^. Charge states are given for two cases, 
one with k evaluated according to the procedure give n above , and o ne with k — 00, i.e. a 
Maxwellian distribution is enforced throughout as in iLamingl (120041 ). Except in cases with 
the lowest values of electron-ion equilibration, a k electron velocity distribution produces less 
ionization than the Maxwellian distribution. With little or no electron-ion equilibration, the 
K, distribution has more electrons above the O"*"^ ionization threshold than the corresponding 
Maxwellian, and so higher ionization states result. At higher electron temperatures, many 
electrons are above the threshold in the Maxwellian, and the k distribution puts more of 
these at even higher energies, where the ionization cross section is starting to decrease. 
Thus at high electron-ion equilibration rates, the k distribution produces less ionization. 
We have previously identified these CcLSGS cLS being necessary to produce the observed degree 
of ionization, and so we conclude that at least in the region of the fast solar wind interior 
to freeze-in, the electrons are heated by lower hybrid waves into a Maxwellian distribution 
function, as would be expected in the case that u >> \/2kvi (see Appendix B). Henceforward, 
K — i> 00 is assumed in all mo dels. The electron power law component observed by Ulysses 
(e.g. iMaksimovic et al.l |2000| ) must be generated by another heating mechanism, possibly 



by lower hybrid waves further out in the solar wind where the ma gnetic field is redu ced so 



that tu « ^/2kvi, or by the damping of whistler turbulence (e.g. IVocks et al.ll2005l ). This 
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conclusion is at variance with the assumption of lEsser fc Edgarl (120001 ) . who argued that 
such a tail on the electron distribution should persist all the way down to the Sun, based on 
th e insensitivity t o such an electron distribution of the elect ron temperatu re diagnostic used 
by lWilhelm et al.l (Il998l ). However, as discussed previously (lLamingll2004l ). other diagnostics 
exist which are sensitive to suprathermal electrons, and these generally reveal insignificant 
populations. Such supr athermal electron tails are often attributed to the Landau damping 
of kinetic Alfven waves (jVinas et al.ll2000l ). which if true, provides no obvious reason for the 
electron heating to be suppressed below 1.5 Rsun- By coupling the electron heating to the 
ion acceleration as in our model, this constraint is automatically satisfied. 



4. Fe Charge States Measured by Ulysses and ACE in Polar Fast Wind and 

the Equatorial Fast Wind 

The creation and destruction of an ion species is governed by the following rate equation 
for its density rii, 

d/Ti ' 

—1 = -Cifune - {Rdr,i + Rrr,i)nine + Ci^iTli^in^ + {Rdr,i+l + -Rr-r,i+l)^i+l^e " 2njf /r (14) 

at 

where rig is the electron density, Cj is the collisional ionization rate, Rdr,i is the dielectric 
recombination rate, RrT,i is the radiative recombination rate, and v and r are the solar 
wind velocity and radius. These the ionization and recombination rates depend on the local 
electron temperature and the speed of the ion. When the electron density is reduced by the 
expansion such that the last term dominates on the right hand side, the charge states are 
said to have frozen- in. Depending on the element concerned, freeze-in occurs between 1.5 
and 3.5 Rq heliocentric distance. Observation of in-situ charge states therefore allows the 
inner corona to be studied in detail. 

We examined charge state distributions of Fe for a variety of fast solar wind periods 
using data we obtained from both ACE and Ulysses. The Ulysses SWIGS data has a time 
resolution of 3 hrs (obtained via private communication, R. von Steiger), while the ACE 
SWIGS data has a time resolution of 1 hr (obtained from the ACE Science Center). The 
coronal holes were identified by examining the solar wind speed, the O^^/O"'"^ ratio, and 
the C"'"^/C"'"^ ratio. The mean charge state distributions were calculated for the 2 — 4 day 
period when either Ulysses or ACE was in a high speed stream not attributed to ICMEs, but 
consistent with coronal hole solar wind. We examine three different types or cases of fast solar 
wind; (1) four periods observed in a single well- formed solar minimum polar coronal hole, 
(2) four periods inside a single solar maximum polar coronal hole, and (3) 14 periods inside 
separate fast solar wind streams from solar maximum equatorial coronal holes observed at 
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ACE in 2005. We will compare the observed charge state distributions with those predicted 
using various heating rates. 

Figure 2 shows Fe charge states observed in the fast solar wind by Ulysses (Case 1). 
There are six time periods, four days in duration, identified as wind originating in the polar 
coronal holes during the southern polar pass. The relative abundance of the charge states 
for these four day periods are averaged and plotted. T he distributions tend to peak between 



Qpe = 10 - 11. This is consistent with the finding of iKo et al.l (119991 ) for the years 1994 



1995. A solution from the model with 'y ' = 0.025 and an initial speed at l.O5i?0 of 10 km s ^ 



(the preferred value from iLamingl |2004| ) is shown to closely match the observed distribution 
function. Fe charge states for a variety of these models are given in Table 2, and 7' is 
defined in section 3. This demonstrates how lower hybrid wave damping can additionally 
heat the solar wind before the charge states freeze-in matching observations, and that only 
one characteristic density scale length is required. 

Figure 3 shows the Fe charge state distributions for a solar maximum coronal hole 
observed by Ulysses SWICS. These four four-day periods occurred during the northern polar 
pass toward the end of 2001. Again the charge states represent the mean charge state 
distribution observed during these periods. The maximum charge state for this period is 
near 10, although a large fraction remains in Qpe = 9. These distributions are compared to 
a model with a heating rate of 7' = 0.015. However it is clear that the observed charge state 
distributions are broader than the model, indicating that a range of density scale lengths 
may be present. 

Figures 4 and 5 show the Fe charge state distribution for fast solar wind associated with 
an equatorial coronal hole observed by ACE SWICS during solar minimum. For these plots, 
14 periods are selected ranging from 2-4 days in duration. These periods represent different 
passes through fast solar wind, whereas the previous two examples represent different periods 
within the same coronal hole. The periods are separated into those that peak at Qi^e = 9 
(Figure 4) and those that peak near Qpe = 10 (Figure 5). The charge state distributions 
that peak at Qp^. = 9 are fit best by a model 7' = 0.01, while the charge state distributions 
that peak at Qfe = 10 are best fit by a model 7' = 0.015, which implies more heating by 
wave damping is needed for these. The model results fail to catch the shoulder of the charge 
state distributions out to a charge state of 14 and beyond. 

Since the ionization model assumes solar minimum conditions over the poles, there 
may be additional characteristics that need to be considered for a full understanding of the 
contribution o f wave damping for t he heating of ions in different solar wind streams from 



coronal holes. iMiralles et al.l (I2OOII ) observed that the O^"*" outflow speed was lower for the 



equatorial coronal holes than for the polar coronal holes, despite the fact that their final solar 
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wind speeds at 1 AU are not significantly different. The wind acceleration must occur further 
out, possibly beyond the radii associated with freeze- in, leading to lower degrees of enhanced 
ionization. It may be that the higher densities in these structures prevent the plasma from 
becoming collisionless further from the sola r surface than is th e case in polar coronal holes, or 
that the increased magnetic field strength (IWang et al.l 120001 ) inhibits the turbulent cascade 
until the field has decreased sufficien tly further out, since three- wave interaction probabilities 
go as {uo/Vtif JLuo k Melroselbood ). 



Ulysses charge states from the solar minimum polar coronal hole during 1994 are ob- 
served to be higher than the charge states in the solar maximum polar coronal hole observed 
by Ulysses in 2001. The charge states in the equatorial coronal hole observed at ACE in 
2005 are even lower, with an average charge state of 9+ for Fe. The trend toward lower 
average charge states correlates with the average of the solar wind speeds during the peri- 
ods measured. The average speed of the solar wind alpha particles observed from the solar 
minimum equatorial coronal hole at ACE was 650 km/s (ranging from 600 - 800 km/s) , 
while solar wind alpha speeds from the solar minimum polar coronal hole at Ulysses were 
at 800 km/s +/- 20 km/s. The wind speeds were more variable in the ACE observations. 
The relationship of solar wind alpha speed and mean Fe charge state and therefore lower 
hybrid damping rate indicates that less additional electron heating is needed in slower fiows 
at lower latitudes. The electron te mperatu r es req uired to produce the observed charge are 
slightly lower than was the case in iLaming) (120041 ) . due to the small downwards revision of 
these charge states. We find the collisionless equilibration parameter 'y' in the range 0.01 - 
0.025, which corresponds to the curves labeled 0.2 - 0.5 in Figure 5 of iLamingl (120041 ). This 
represents a small chang e , bec oming a little more consistent with the observational error 
bars from IWilhelm et al.l (119981 ). 



5. Discussion: Relation to MHD Turbulence 

The cascade of Alfvenic MHD turbulence from small to large is expected to pro- 
duce cross field v ariati ons on small length scales of density, magnetic field and velocity. 



Markovskii et al.l (120061 ) consider in most detail the velocity shear produced by perpendic- 
ularly propagating Alfvenic fiuctuations. Such perturbations will also produce small scale 
perpendicular variations in magnetic field. While magnetic curvature may also give rise to 
lower hybrid waves, it is the longitudinal variation that is important in this respect. However 
k\\ does not cascade to large values, unlike so we are left with cross-field density gradients 
as the most effective means of generating lower hybrid waves. Significant density fiuctuations 
cannot arise solely from Alfvenic turbulence. Waves of a compressive character such as fast 
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mode waves are required. iLuo fc Melrosd (120061 ) . in an analytic consideration of three- wave 
processes involving two or three Alfven waves and at most one fast mode wave, show that a 
critical frequency, cUc, exists, such that below this frequency fast mode waves dominate, and 
above Alfven waves are more prevalent. In the simplest case involving one fast mode wave 
and two Alfven waves, with one of the Alfven wave frequencies, uji, much lower than the other 
two, Uc ~ 2ujp(^ujl^^ , putting Wp ~ Wa in their equation 25. For a density of 10^ cm~^ and 
magnetic field of 1 G, cj; ~ 0.07 which corresponds to a wave period of 90 s, assuming that 
in our case, fast mode turbulence must dominate up to frequencies correspondi ng to the ion 



gyrofrequency. Waves of t his period or longer are reported in coronal holes, (jOfman et al. 



19991 : iBanerjee et al.ll2001), and so the existence of fast mode turbulence would not seem to 
be a problem. IChandraru (120051 ). adopting a numerical approach to consider a wider range of 
processes involving Alfven and fast n aode waves, finds different regions of fc- space dominated 
by Alfven or fast mode turbulence. ISuzuki. Lazarian. fc BeresnyakI (120061 ) conjecture that 
the degree of couphng bet ween fast mode s and Alfve n modes should b e reduced in strong 
Alfvenic turbulence; both iLuo fc Melrosd (120061 ) and IChandranI (120051 ) only consider weak 
Alfvenic tur bulence. Slow mode waves have not been considered as part of the turbulent 
cascade, but iBallai et al.l (120031 ) discuss the behavior of slow mode solitary waves, which 
could also produce suitable density gradients in perpendicular propagation. 

Direct evidence of turbulent cascade down to these small scales has been hard to come 



by. However Cluster (IBale et al.ll2005l ) has observed electric and magnetic fluctuation spec- 
tra in the solar wind near the terrestrial bowshock. In the inertial range, both spectra 
resemble a Kolmogorov k~^^^ behavior. However as the dissipation range is approached, at 
kvg ~ 1 where Vg is the ion gyroradius, the magnetic spectrum steepens, while the electric 
spectrum flattens. Clearly the turbulence is becoming more electrostatic. Measurements of 
the phase speed suggest the turbulence is mainly Alfvenic in nature, transforming to kinetic 
Alfven waves in the dissipation range and damping by heating the ambient plasma as the 
waves become more electrostatic. This is entirely consistent with the approach taken in this 
paper, except that the electrostatic waves we consider derive not directly from the cascading 
turbulence, but ClS db by-product of the damping. 

We determine a likely order of magnitude of 7' ~ 0.01 from charge states observed in 
the fast solar wind. Assuming a wave electric field given by equation A3 and a Maxwellian 
distribution for the a parti cles that excite the waves, we determine a density scale length 
L ~ Avg from figure 11 in iLamingl (120041 ). where Vg is the gyroradius of the a particles. 
Observations of a partic les in the solar wind typically give a k distribution in the parallel 
direction with k ~ 3 — 6 (jCoUier et al.lll996l : IChotoo et al.lll998l ). which if it has an origin in 
energization by ion cyclotron waves would correspond to k ~ 3.5 — 6.5 in the perpendicular 
direction. This would increase our density scale length estimate to L ~ 6 — 10r„. We showed 
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earlier that ion heating is neghgible compared to electron heating for lower hybrid waves 
excited under solar wind parameters. Thus the electric field estimate in equation A3 may 
be exceeded by as much as an order of magnitude. Hence 7' ~ 10"'*, and the density scale 
length becomes L ~ 10 — 25rg. We have no direct knowledge of the a particle gyroradius. 
The perpendicular velocities for protons and O VI are ~ 200 and ~ 300 — 400 km s~^ at a 
heliocentric distance of about 2 Rq, where the magnetic field is ~ 1 G and the particle density 
is ~ 10^ cm~^. The proton gyroradius is then 0.02 km, and the proton inertial length is 0.2 
km. At 1 AU, the minor ions fiow about an Alfven speed faster than the protons, presumably 
refiecting enhanced perpendicular heating b y ion cyclotron waves c loser to the sun. The a 



particles also fiow faste r than the protons (IReisenfeld et al.l I2OOII : Ivon Steiger et al.l 12000 



Neugebauer et al.l 119961 ). so we will assume an a particle perpendicular velocity of 400 km 



s~^, similar to O VI. This gives an a particle gyroradius of 0.08 km, and the range of density 
scale lengths corresponds to 0.8 - 2 km, or 4 - 1 proton inertial lengths . This is very similar 
to the scale length of velocity shear discussed by lMarkovskii et al.l (120061 ). Thus it seems that 
our results support ideas of turbulent cascade of MHD fiuctuations in the solar wind. Other 
observat i ons o f density structures observed by Interplanetary Scintillation are reviewed in 
Laminel kooi \. 



Another often cited mechanism for electron heating in the solar ^y ind is the Landau 



damping of kinetic Alfven waves (e.g. ICranmer &: van Ballegooijenll2003l ) by the small par- 



allel component of the wave elect ric field. Th i s proc ess, and that for fast mode waves also, 
are considered in some detail by IVinas et al.l (I2OOOI ). with application to the formation of 
suprathermal electrons in the chromosphere. Their solution of the dispersion relation only 
yields significant parallel electric fields for wavevectors approaching 10c/ Upi for both fast 
mode and kinetic Alfven waves. At ~ c/upi, the parallel electric field is a factor ~ 10~^ 
smaller than the perpendicular field, and smaller by at least an order of magnitude than the 
parallel electric field due to lower hybrid waves if the density gradients we infer above are 
related to Alfven or fast mode w avevectors. Mor e generally, one can estimate the parallel 
electric field from the expression (lHasegawalll976l . equation 37) 



Ell 



(15) 



which would suggest that wavelen gths of order the pr oton gyroradius are required. This is 
less restrictive than the results of IVinas et al.l (120001 ). but still significantly more so than 



what we determine for lower hybrid waves. It would appear then that electron heating 
by kinetic Alfven waves places even more stringent requirements on the turbulent cascade 
than do the lower hybrid waves. Another potential problem is that no natural means of 
suppressing electron heating at altitudes within 1.5 Rq heliocentric distanc e exists, which 



might lead to confiict with the SUMER electron temperature observations of IWilhelm et al. 
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(119981) and others. 



Conclusions 



We have revisited the model of iLaming) (12004J ) for increased ionization in the fast solar 
wind with the analysis of new data, and some refinements to the theory. The important new 
development, due in part to the revision downwards of ch arge states obs erved in the solar 
wind, and more significantly to the correction of an error in iLamina (120041 ) in the expression 
for the wave energy density, is that the degree of density inhomogeneity required in the 
fast solar wind acceleration region to produce the required heating is now reduced from a 
length on the order of an ion gyroradius, to be in the range 10 - 25 gyroradii. As discussed 
above, this makes a connection between these density gradients and the outcome of the 
MHD turbulent cascade much more likely. Our finding that ion excited lower hybrid waves 
develop stronger parallel ele c tric fi elds than kinetic Alfven waves reinforces suggestions like 
that of iBegelman fc Chiuehl (119881 ) that such mechanisms may have more general relevance 
elsewhere in astrophysics. 

In our models, the appropriate choice of the electron-ion equilibration parameter can 
yield a fit for the observed distribution functions for various fast solar wind streams. This 
supports the idea that existence of density gradients from MHD turbulent cascades can allow 
lower hybrid waves to oscillate and damp, imparting their energy to electron heating. This 
electron heating due to the damping of the lower hybrid waves occurs in the region beyond 
1.5i?0, where ions begin to execute large gyro- orbits under the influence of perpendicular 
heating by ion cyclotron waves. The solar wind flows with the highest alpha particle speeds 
also have the highest average charge states, implying that the electron heating is most 
efficient in fast outflow from polar coronal holes during solar minimum. Lower fast wind 
speeds result in lower average charge states of Fe, clearly indicating a connection between 
electron heating and ion acceleration. Finally, in Table 3, we give predicted charge states 
for Ne, Mg, Si, S, Ar, and Ca, again for an initial wind speed of 10 km s^^ at l.O5i?0. 



This work has been supported by NASA LWS Grant NNH05AA05I (JML) and by an 
NSF SHINE Postdoctoral Fellowship ATM-0523998 (STL). 
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A. Wave Trapping and Threshold Electric Field 

The equation of motion for a particle in the electric field of a wave is 

dv qE qE 

— = — exp tip = — cos ip ( Al j 

at m m 

where ip = uot—kz and real parts are taken in the last step. Writing d'^z/dt^ = — (l/fc) d'^ip/dt^ 
qE / m cos ip we arrive at 

92 di> qEkdip 2 9ip . , . 

d^m = ^ 9^ ^^"^ = -T^sm^ (A2) 

where ujt = \J —qEk/m is the particle trapping frequency. A particle is said to be "trapped" 
when V > {uj ~ ujt) /k during its cyclotron orbit. The time spent in the trapping region 
during a cyclotron orbit is 26/Q where cos^^ = 1 — ujt/^ ~ 1 — 9'^ /2, so the time is 2/VL x 
^JtLjtJu)- The trapping is defined as effective when the time spent in the trapping region is 
> l/cuy, so 8ti;|. = fi^cij, which then gives 



Bo 4 I y kc 



(A3) 



which is the result of iKarneyi (1l978l ). derived by integration of the Hamiltonian equations of 
motion. 



B. Electron Distributions from Lower-Hybrid Wave Heating 
We write the Boltzmann equation as 



df 1 d 



dt v± dv± 



OV\\ OV \ OV\\ ' 



(Bl) 

where A_[_ and A\\ are coefficients of dynamical friction and -D^_|_, D\\\\, -^J-lh ^"^^ 
velocity diffusion coefficients. Integrating the Boltzmann equation over v±_ from to 00, we 
have 

^ j f2Txv^dvx. + j Di\±-^27rv±dv± + D\\\\-^ j f2nv±dv±'^ = (B2) 
in steady state conditions (i.e. df /dt = 0). This integrates to 

j f2Tivi_dvi_ cx exp j A\\/D\\\\dv^ = exp {-vl/2v^\\) (B3) 
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when A\\ = fy/r, D\\ 
the stopping time, 



v^^/ts, with the deflection time Td » Ts 



= Vt\\Vt±/Td and D\\ 

27mef|/3a;pgln A for fast electrons in cold plasma, with In A ~ 30 
being the Coulomb logarithm. For comparison, the electron-electron equilibration time. 



"eg 



A/2/37r^/^?T,ef ^/cUpg In A, is a factor of a few higher. 

For electron heating by L andau damping , the non-zero diffusion coefficient and coeffi- 
cient of dynamical friction are (jMelroselll986l . equation 10.83 and 10.90 with s = 0) 

(B4) 



All 



SttV i? 



e- V 
e- V 



where in the small gyroradius limit V = (O, 0, v\\) and e is the wave polarization vector. The 
number density of wave quanta is given by A^, R is the ratio of electric energy to total energy 
in the wave, such that R J Nhwd^k/ (27r)'^ = 6E'^/87t. For lower hybrid waves 



1 + T5r sm^ 9 ^ + 

UJ'' kvf 



(B5) 



where (x) is the usual plasma dispersion function and Vi = ^J'^oT^Jmi is the ion thermal 
velocity. We consider two possibilities. When uj » ^/2kvi, 0^1 — k'^vf/uj'^ and 



UJ 



+ cos^ e 



and when uj « ^/2kvi, (j) ~ uj"^ /k'^vf^ and 



UJ 



ujp^ cos^ 9 



l + ulsinH/ni + ujlJk^vr 



(B6) 



(B7) 



In either case, the ratio of electric to total energy in the wave R = {ujd^K^/duj) ~ nl/2uj:^ 



and 



e- V 



= f y cos^ 6'. For UJ « \^kvi, both uj and the growth rate 7 vary as cos^ 6*. 
Consequently wave generation will occur mainly at cos 6* = ujpi/ujpe, beyond which lower 
hybrid waves are progressively Landau damped by parallel propagating electrons. This then 
yields 



All 



9fc|| \ me / "^pe""]! 



tel. 



:b8) 



where d/dk\\ {q6E/me)'^ is evaluated at k\\ = uj /v\\. For all cases of interest, A\\ « v\\/t and 



so 



/ ~ exp j v\\/ [v\ + rD||||] dv\^ = [v% + vy2K] " 



(B9) 
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where K = \4:'Kd/dk\\ {q5E/mef x vrtUpjrie/Sco'pg In A] . We have assumed that 9/(9fc|| {q6E/ 
is independent of fcy at uj/v\\ and hence that D\\\\ oc l/fy in the final step. 

Where u; >> \^kvi, both lu and 7 remain non-zero as cos 6' ^ 0, so wave generation 
occurs for all cos 6' in the range —Upi/upe Upi/upe. In this case then 



,2 

rUe 



2 2 



III ofe|| y me J K^tiji 

All = ^ jT' , (BIO) 

II teX 

and the electron distribution function integrates to a Maxwellian 

/ {vii) oc exp (^-vl/2/ + {7^71^00^6^^^ InAP) And/dkn {q6E/ 
= exp (^-vl/2/ + (ujyevjj {rs-r^/v\i) d/dk\\ {q6E/m,) 

We also estimate Aii ^ lO^'^^/ne << vh/ts as required with B = lOG. 



,2 



C. Atomic Data Updates 



We have also updated the atomic data for dielectronic recombination of K-shell, and 
L-shell ions. Recomb i nation from H- to He-like and from He- to Li-like are taken from 
Dasgupta fc Whitneyl (120041 ). T he successive isoelectronic sequence s Li-, Be-, B-, C-, N - 
, 0-, and F -l ike ar e taken from Colgan. Pindzqla. Sz Badnell (|2004 ). Colgan et al. ( 20031) . 



Alturi et all fl2004ri. IZatsarinnv et all (l2004al ). iMitnik fc Badnelll (120041 ). IZatsarinnv et al. 
( 2003! ). and Gu ( 20031 ) respectively. Additional ly dielelectronic r e combin atio n fro i n Ne- to 



Na-like and from Na- to Mg-like are taken from IZatsarinny et al.l (j2004bl ) and iGul (120041 ). 
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Table 1: 0+6/Q+7 Abundance Ratios 
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Table 2: Fe Charge State Fractions 
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Table 3: Minor Ion Charge State Fractions, Vstart = 10 km s ^ 
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Fig. 1. — Approximate (narrow solid line, from equation 12) and true (wide solid line) k 
distributions for varying k. The underlying Maxwellian is shown in each case as a dashed 
line. 
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Coronal Holes in 1994 Measured by Ulysses SWICS 
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Fig. 2. — Observed (colored solid lines) and modeled (dashed line) Fe charge states in fast 
wind from polar coronal holes observed at solar minimum by Ulysses SWICS in 1994. 
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Coronal Holes in 2001 Measured by Ulysses SWICS 
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Fig. 3. — Observed (colored solid lines) and modeled (dashed line) Fe charge states in fast 
wind from polar coronal holes observed at solar maximum by Ulysses SWICS in 2001. 
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Coronal Holes in 2005 Measured by ACE SWICS 
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Fig. 4. — Observed (colored solid lines) and modeled (dashed line) Fe charge states in fast 
wind peaking at Q^e = 9 from equatorial coronal holes observed by ACE in 2005. The average 
charge states are lower still than the solar maximum coronal holes, and the distributions are 
generally broader. The a particle wind speed is also the lowest in these cases. 
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Coronal Holes in 2005 Measured by ACE SWICS 
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Fig. 5. — Observed (colored solid lines) and modeled (dashed line) Fe charge states in 
fast wind peaking at Qpe = 10 from equatorial coronal holes observed by ACE in 2005. 
The average charge states are lower still than the solar maximum coronal holes, and the 
distributions are generally broader. The a particle wind speed is also the lowest in these 



cases. 



